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Abstract 

We study the statistics of the recurrence times between earthquakes above a cer- 
tain magnitude M in California. We find that the distribution of the recurrence 
times strongly depends on the previous recurrence time tq. As a consequence, the 
conditional mean recurrence time t(tq) between two events increases monotonically 
with To- For To well below the average recurrence time f, t(tq) is smaller than T, 
while for To > T, t(to) is greater than r. Also the mean residual time until the next 
earthquake does not depend only on the elapsed time, but also strongly on to- The 
larger To is, the larger is the mean residual time. The above features should be taken 
into account in any earthquake prognosis. 



Recently, Corral [1,2] studied the recurrence of earthquakes above a certain 
magnitude threshold M in spatial areas delimited by a window of L degrees 
in longitude and L degrees in latitude. He found that the distribution D{r) of 
recurrence times r scales with the mean recurrence time r as 

D{r) = \f{r/r), 

T 

where the function /(B) is quite universal and independent of M. For B below 
1, / can be approximated by a power-law, while for B ^> 1, /(B) decays 
exponentially with B. As a consequence of the deviation from a Poissonian 
decay [3], the mean residual time to the next event increases with the elapsed 
time [1,2]. 

In this paper, we study the statistics of the recurrence intervals of the Cali- 
fornia database [4] and find that both quantities, the recurrence interval dis- 
tribution D(t) and the mean residual time to the next earthquake strongly 
depend on the previous recurrence time interval t . 
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We study the records from the local Californian earthquake catalog [4] for the 
period 1981-2003 in the area 30.5-38.5N latitude, 114-122W longitude, with 
minimal magnitude threshold value M = 2 and minimal recurrence times 2 
mins [5]. Similarly to Corral [1], we consider the earthquakes in the region 
that are above a certain threshold M, as a linear process in time {£;} without 
taking into account the spatial coordinates of the event hypocenters. We are 
interested in the recurrence intervals Tj = ti — tj_i between these earthquakes. 

In records without memory, the (conditional) distribution D(t\t ) of recur- 
rence intervals r that directly follow a certain interval To, does not depend on 
the value of To and is identical to D(t). In contrast, in records with long-term 
memory, there is a pronounced dependence of D(t\tq) on r [7,8]. To study 
possible memory effects in the earthquake records with a reliable statistics, we 
have studied the conditional distribution D(t\tq) not for a specific r value, 
but for values of r in certain intervals. To this end, we have sorted the record 
of N recurrence intervals in increasing order and divided it into four subrecords 
Qi, Q2, Q3 and Q4, such that each subrecord contains one quarter of the to- 
tal number of recurrence intervals. By definition, the N/4 lowest recurrence 
intervals are in Qi, while the N/4 largest intervals are in Q4. 

Figure 1 shows D(t\tq) for To averaged over Qi and Q±. For comparison, we also 
show the unconditional distribution function D(t). To improve the statistics, 
we used logarithmic binning. We considered time scales from 2 minutes to 
10t, with 50 log-bins, counted the number of recurrence intervals within each 
bin and divided it by the size of the bin. To further improve the statistics, we 
averaged the probability distribution over threshold values M = 2.25 ... 2.75 
around M ~ 2.5. Finally, we normalized the probability distribution to obtain 
the probability densities of interest. The figure shows that for r well below 
its mean value r, the probability of finding r below (above) r is enhanced 
(decreased) compared with D(t) for tq in Qi, while the opposite occurs for tq 
in Q 4 . 

By definition, r(r ) is the mean recurrence intervals, when the two events 
before were separated by an interval To- The memory effect in the conditional 
distribution function D(t\t ) leads to an explicit dependence of f (r ) on r . To 
calculate f (to), we divided the sorted (in increasing order) record of recurrence 
intervals into 8 consecutive octaves. Each octave contains iV/8 intervals. In 
each interval, we calculate the mean value. We studied f as a function of t /t, 
where now To denotes the mean recurrence time in the octave. Figure 2 shows 
t(tq)/t as a function of tq/t and clearly demonstrates the strong effect of the 
memory. Small and large recurrence intervals are more likely to be followed 
by small and large ones, respectively, f/r is well below (above) one for To/t 
well below (above) one. When the recurrence intervals are randomly shuffled 
(no memory), we obtain f(r )/r = 1, see Fig. 2, open symbols. 
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A more general quantity is the expected residual time f(x|ro) to the next 
event, when time x has been already elapsed. For x = 0, f(0|r ) is identical 
to t(t ). In general, t(x\t ) is related to D(t\t q ) by 



oo ,00 

f(x|r ) = J(t - x)D{t\t ) dr I J D(r\r ) dr 



For uncorrelated records, D(t\t ) is Poissonian, and t(x\to)/t = 1. 

Figure 3 clearly shows that t(x\tq) depends on both x and To. With increas- 
ing x, the expected residual time to the next event increases, as is shown in 
Fig. 3, for values of r from Qi and Q A (top and bottom curves). Thus, when 
t increases, t(x\t ) increases for all values of x. The middle curve shows the 
expected residual time averaged over all To, i.e. the unconditional residual time 
f(x). In this case, the interval between the last two events is not taken into ac- 
count, and the slower-than-Poisson-decrease of the unconditional distribution 
function D(t), Eq. (1), leads to the anomalous increase of the mean residual 
time with the elapsed time [3]. 

Our results for the unconditional residual time function for Californian earth- 
quakes are very similar to the results of Corral [6] that were obtained for 
worldwide earthquake records. As shown here, there is a strong memory in 
the earthquake recurrence intervals, which influences significantly the resid- 
ual time. Similar memory effects have been obtained recently for river flux, 
temperature and precipitation records (see [7,8]). 

To summarize, we have studied the memory effect in the earthquake events and 
showed that the distribution of the recurrence times and the mean residual 
time until the next earthquake strongly depend on the previous recurrence 
time. The conditional mean recurrence time between two events monotonically 
increases with tq. These results should be taken into account in an efficient risk 
evaluation and forecasting of earthquakes. It is very plausible that the origin 
of these effects is due to long-term persistence in the earthquake occurence. 
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Fig. 1. Conditional probability distribution for the recurrence time intervals between 
earthquakes above a threshold M ~ 2.5 following recurrence time To from the first 
quarter (circles) and the last quarter (squares) of the recurrence time, and the 
unconditional probability (stars). To improve statistics, averages were taken for 
2.25 < M < 2.75. 
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Fig. 2. Expected recurrence time t(tq) between earthquakes above thresholds 
M ~ 2.5 (full circles), M ~ 3 (full squares), M ~ 3.5 (full triangles up), and 
M ~ 4 (full triangle down) following To taken from the eight octaves described 
in the text. Averages are taken in intervals M ± 0.25 to obtain better statistics. 
The open symbols represent the analysis of the randomly shuffled recurrence time 
record, yielding f (tq)/t ~ 1. 
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Fig. 3. Conditional mean residual time to the next earthquake above a thresh- 
old M ~ 2.5 following recurrence time To taken from the first quarter (bottom 
curve) and the last quarter (top curve) of the recurrence intervals, and uncondi- 
tional mean residual time (middle curve). To improve statistics, average is taken for 
2.25 < M < 2.75. 
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